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General phenomenological theory of hydrodynamic waves in regions with smooth loss of convexity of isen- 
tropes is developed based on the fact that for most media these regions in p -V plane are anomalously small. 
Accordingly the waves are usually weak and can be described in the manner analogous to that for weak shock 
waves of compression. The corresponding generalized Burgers equation is derived and analyzed. The exact 
solution of the equation for steady shock waves of rarefaction is obtained and discusses. 



PACS numbers: 47.40.-x, 43.25.+y, 52.35.Mw 

Extremely great importance of the so-called elementary 
wave patterns, such as shock and rarefaction waves, in hy- 
drodynamics has been known since the pioneering studies of 
Riemann [[]]] carried out in the beginning of 1860s. The im- 
portance is generally connected with (i) possibility to picture 
general fluids flows as nonlinear superposition of elementary 
waves, propagating as separate entities and (ii) reflection with 
the elementary waves the asymptotic behavior of general so- 
lutions of the set of hydrodynamic equations. For this reason a 
great deal of attention has been payed to study the elementary 
waves, see, e.g., the review article |Q| and references therein. 
Nonetheless, a certain type of the elementary waves, namely 
shock waves of rarefaction and sonic waves [0], which admits 
general analytical consideration, has not been considered in 
such a manner. The gap is filled up in the present Letter. 

It is well known that the sign of variation of pressure and 
density in a weak shock wave is associated with the sign of 
(d 2 V/dp 2 ) s . Usually the sign of the derivative is positive, 
and accordingly the shock wave is compressive. However, in 
some specific areas of p -V plane, e.g., close to a critical point 
but still in a single-phase region, (d 2 V/dp 2 ) s may become 
negative. The phenomenon is known as smooth loss of con- 
vexity of isentropes [|J] and is quite general for many types of 
equations of states of real materials. Even the van der Waals 
equation possesses a region around the vapor dome, where 
(d 2 V/dp 2 ) s < 0. In 1942 Bethe showed that all fluids 



having sufficiently large specific heats exhibit smooth loss of 
convexity of isentropes in the neighborhood of the saturated 
vapor line [Q] . 

The smooth loss of convexity of isentropes creates pre- 
requisites for a shock wave of rarefaction (SWR) to come 
into being. Moreover, other types of anomalous wave struc- 
tures, such as sonic shocks, whose Mach number(s) before 
the shock, after it, or both is (are) identically equal 1, may 
be observed in this region and its vicinity, see Refs. [|z|3|] for 
more details. Bearing in mind all mentioned above, it is highly 
desirable to have a certain universal equation, which is not 
connected with any particular equation of state but neverthe- 



less can describe hydrodynamics waves in the discussed re- 
gion. Despite the apparent importance of the problem and its 
long history, much to our surprise, nobody has payed attention 
to the following consequence of the fact that usually the re- 
gion with the smooth loss of convexity of isentropes occupies 
a small domain in p-V plane. The consequence is that the 
corresponding waves are always weak (otherwise the wave 
curves [||] go far away from the domain), and therefore the 
problem of the wave description is tractable within the frame- 
work of the macroscopic hydrodynamics in the manner anal- 
ogous to that for weak shock waves of compression. To be 
precise we must stipulate that we do not consider any singu- 
larities related to the very vicinity to the critical point. In other 
words, it is implied the system is far enough from it to avoid 
any effect of the singularities. 

In the present Letter the mentioned approach is applied to 
the problem. As a result the generalized Burgers equation, 
which describes hydrodynamic waves in the region, where 
the derivative (d 2 V/dp 2 ) s is small and may change sign is 
derived. The equation is valid for both signs of the derivative 
and therefore describes all the variety of waves, which may 
exist in this region. As an example the equation is employed 
to study steady shock waves of rarefaction (SWR). The gen- 
eral exact solutions of the equation describing such shocks is 
obtained and analysed. 

To begin with, let us consider the conventional Burgers 
equation, which governs a weak shock wave of compression. 
In the laboratory coordinate frame for a wave advancing in the 
negative direction of x-axis the equation reads [0] 
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where p' is the pressure variation, so that the pressure in the 
wave profile p equals to p\ + p' (here and in what follows 
subscript 1 indicates the value of the corresponding quantity 
in the unperturbed medium, i.e., at x — > — oo, while subscript 
2 will stand for the value of the same quantity behind a wave, 
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at x — > oo), c(p) is the velocity of sound, 
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2V 2 V dp 2 



V denotes the specific volume (V = 1/p), a is a dissipative 
constant [pi] 
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?7 and £ stand for the coefficients of viscosity, k for the thermal 
conductivity and c V:P for the corresponding specific heats. 

Eq. ([!]) may be regarded as leading approximation to expan- 
sion of a certain more general nonlinear dissipative equation 
in powers of both weak nonlinearity and weak dissipation. It 
is important that the nonlinear term in this approximation is 
considered in non-dissipative limit. Dissipative corrections to 
it have additional smallness and may be neglected. 

In our case | (d 2 V/dp 2 ) s | is anomalously small (it even can 
vanish). Accordingly, the lowest nonlinearity ap' (dp' / dx) 
may not correspond to the leading nonlinear term any more. 
It means that higher order (in p') terms must be taken into ac- 
count. On the other hand, the region, where (d 2 V / dp 2 ) s < 0, 
is narrow, therefore the third derivative (d 3 V / dp 3 ) s also has 
a certain smallness in this region. Only the fourth derivative 
(d A V/dp A ) s does not have any smallness there. Thus, at small 
but finite p' terms O ((p 1 ) 2 ) , O ({p') 3 ) and O ((p') i ) may 
yield contributions of the same order of magnitude and must 
be taken into account simultaneously. It is important that dis- 
sipative corrections to the nonlinear terms remain negligible 
and all such terms may be considered in the non-dissipative 
limit, in the same manner as it is for the conventional Burgers 
equation. All the above mentioned is related exclusively to 
nonlinear terms. Regarding the dissipative term on the right 
hand side of Eq. (|l]), there are no specific reasons for it to be 
small and no correction to this term is required. 

To derive the desired nonlinear corrections to the left hand 
side ofEq. ^ let us consider an arbitrary adiabatic flow in the 
form of a traveling wave. Due to the adiabatic condition the 
profile of entropy of such a flow should not change, and since 
the state of the medium ahead of the shock front is spatially 
homogeneous with s — const, we obtain that s = const 
anywhere. Then, we may employ the general solution of the 
continuity and Euler equations valid for simple waves [^j and 
reduce these equations to the following single equation 



where the flow velocity v is given by the expression 
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v = ± 



dp 



c(p)p(p) ' 
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Sign plus in Eqs. ([z|)-(|3|) corresponds to the wave advancing 
in the positive x direction, minus to negative. According to 



our choice of the direction of the wave propagation in what 
follows we take sign minus. 

It is important to stress that under the above-specified con- 
ditions the reduction of the set of hydrodynamic equations to 
Eqs. (^])-([3]) is an exact result valid for any nonlinear depen- 
dence c(p) and p(p), which may be obtained for the isentropic 
wave. However, if p 1 is small one can expand these functions 
about the point p' = 0. It is seen straightforwardly that if one 
truncates the expansion at zero term, it reduces Eqs. (^)-(Q) 
to the linear acoustic equation. The truncation at the term of 
order 0(p') yields the left hand side of Eq. ([!]). Increasing the 
order of the truncation one can obtain the desired corrections 
to the conventional Burgers equation with any given accuracy. 

Thus, the generalized Burgers equation must have the form 



dp' ... dp' „ d 2 p' 
_ +M(p) _ =Cia _ 
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where u stands for the velocity of a point in the wave profile. 
For our choice of the propagation direction u(p') = v(p') — 

c(p')i- 

To get the equation in the explicit form one must expand 
u(p') in powers of small p' retaining the terms up to O ((p') 3 ) ■ 
However, before doing that it is worth analyzing some general 
features of Eq. (Q). 

First of all let us obtain the general expression for the ve- 
locity v s of a steady traveling shock, when p'(x, t) — > p'(x + 
v s t) = p'(£). In this case Eq. (Q) may be easily integrated. 
Then, taking into account the boundary conditions (p' — > 
0, dp'/d£ -> at f -> -00 and p' -> p 2 - pi, dp'/d£ -► 
at £ — * 00), we arrive at the following expression for v s 
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where (. . .) denotes average overp' [we remind that u = v — 
c < see also Eq. (g)]. 

Let us show now that v s given by Eq. for SWR satisfies 
the evolutionary conditions. The conditions say that in the co- 
moving coordinate frame the velocity of the flow before the 
shock should be bigger than c\, while the velocity behind the 
shock should be smaller than c 2 , i.e., M\ > 1 > Af 2 > 0, 
where M\ 2 stand for the corresponding Mach numbers. In 
the laboratory coordinate frame the medium before the shock 
is in the rest state with zero velocity. Accordingly, in the co- 
moving coordinate frame the flow velocity here is v s . Thus, 
the first evolutionary condition says v s > c\. Behind the 
shock the flow velocity in the laboratory frame is V2- Con- 
sequently, the second condition yields v s + v 2 < c 2 . Bearing 
in mind Eq. (|J) both the conditions may be written as follows 



(-1 



-Ui < -(u) < c 2 - v 2 



-u 2 



(6) 



Finally, taking into account that for a simple wave and the 
chosen direction of propagation the sign of derivative du/dp' 
is opposite to that for (d 2 V/dp 2 ) s [H], i.e., for negative 
(d 2 V/dp 2 ) s function —u(p') monotonically increases with 
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decrease of p', and that for SWR p 2 is smaller than p\, we 
reduce inequality (||) to evident. 

Let us proceed with the derivation of the generalized Burg- 
ers equation. From all the above-mentioned it is clear that in 
the discussed region the derivative (d 2 V/dp 2 ) s may be ap- 
proximated as follows 
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Here p m is the value of p, corresponding to the local min- 
imum of the derivative, e is a small dimensionless quantity 
/i = O(l) and power u 6 as well as the numerical coefficient 
are introduced for convenience of further notations. 

According to Eq (Q), (d 2 V/dp 2 ) s < at p m - A < p < 
p m + A, where A = 2pic 2 e/ fj, 3 = O(e). Then, expanding 
u{p') in powers of p', taking into account that for the prob- 
lem in question p' is of order e, or smaller and dropping terms 
higher than 0(e 3 ), after some calculations we reduce the gen- 
eral equation (|j) to the form of Eq. (jl|), where term ap' should 
be replaced by the following expression 
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and all derivatives (d n V/ dp n ) s are taken at p = pi . 

It is convenient to rewrite the equation in more universal 
dimensionless form. Let us introduce new variables 



P' 



t; C 



(x + Cit). 



In these variables involving expression (0) one can reduce 
the generalized Burgers equation to the following final form 
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where 



f(y) = (z 2 -l)y + zy 2 + V-; z 



Pi - Pr, 



(8) 



(9) 



In what follows we examine solutions of this equation, 
which correspond to SWR. In this case we should supple- 
ment Eq. (|j) with the boundary conditions dy/dC, — > at 
C, — > ±oo; 2/ — *■ 2/1 = at £ — > — oo; y — > y 2 — const < 
at £ — > oo and bear in mind that y < 0; — 1 < z < 1. To 
study steady SWR it is convenient to introduce a traveling co- 
ordinate r/ = ( + vt, where in accordance with Eq (^) the 
dimensionless velocity v equals the following expression 

2 y 3 12 
Then, integration of Eq. (ft) yields 
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2/2(1/3 - y2)(y4 - 2/2) 

ln|2/ - 2/4 1 



2/3(2/3 - 2/4) (2/3 - 2/2) 2/4(2/4 - 2/3) (2/4 

where 2/3,4 are the roots of the equation 
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(10) 



2/ 2 + (4* + 2/2)2/ + 2/ 2 2 + 4^2/2 + 6(z 2 - 1) = (11) 

Let us show that 2/3,4 are always real and one of these root is 
smaller than y 2 , while the other is bigger than 2/1 (we remind 
that by definition 2/1 = 0). In other words the left hand side 
of Eq. is negative at 2/2 < 2/ < and any possible values 
of z and 2/2 ■ To prove it note that for SWR pressure p 2 must 
satisfy obvious conditions p m — A < p 2 < Pi, which in 
the dimensionless variables are transformed into the following 
inequalities 

+ < 2/2 < 0. 



Due to the fact that y 2 enters into Eq. (IT) with a positive 
coefficient to prove negativeness of this polynomial at 2/2 < 
y < it suffices to examine its values at the marginal points 

y = 2/1,2- At 2/ = J/1 = we obtain 

.2 



y 2 +4zy 2 + 6(z' z -l). 



(12) 



For the same reason it is sufficient to inspect the values of 
polynomial (jl/J) at the marginal values of 2/2, namely at 2/2 = 
— (z+1), 2/2 = 0. It is straightforward to see that the marginal 
values of E q. ( Oh are negative at \z\ < 1. Negativeness of 
polynomial ([11 ) at y = y 2 is proved in the same manner. The 
proved relative position of points 2/1.2,3.4 guarantees that the 
derivative dy/dn is negative at 2/2 < 2/ < 2/i f° r anv possible 
values of z and 2/2, i e., the profile of the steady SWR is a 
monotonically decreasing function of rj. 

To end this Letter we present several particular versions of 
general solution (10), when the dependence 2/(77) may be ob- 
tained explicitly. 

(i) z = 1 2/2 = — 2, which corresponds the maximal possi- 
ble amplitude of SWR (pi = p m + A, p 2 — p m — A). 
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(ii) z = 1 1 2/2 1 < 1 Pi = Pm + A, pi - p 2 < A. In this 
case the leading approximation in small \y 2 \ yields 



y = 2/2 



1 + tanh(2/|?y/3) 



(iii) z = -(1 + 2/2), I J/2 1 < 1 [Pi = Pm - A - (pi - 
P2), Pi — P2 *C A]. The leading approximation in I2/2I in this 
case results in the following profile 



V = 2/2 
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1 - tanh (y 2 r)/3) 



Note, that while in the present Letter only steady solutions 
of Eq. (||)-(||) corresponding to SWR are discussed, the equa- 
tion itself describes much more broad spectrum of problems 
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of steady and non-steady flows, including such very interest- 
ing issues, as sonic waves, evolution of arbitrary initial pro- 
files, collision of shocks, etc. The corresponding study is in 
progress and results will be reported elsewhere. 
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